1 



RIEMANN EXTENSIONS IN THEORY OF THE FIRST 
ORDER SYSTEMS OF DIFFERENTIAL EQUATIONS 

^ ■ Valery Dryuma 

O 



X 



Institute of Mathematics and Informatics AS Moldova, Kishinev 
E-mail: valery@dryuma.com; cainar@mail.md 
Abstract 



•*-> 
O 

O 

in 

(N 

O 

. The properties of the Riemann extensions of nonriemannian spaces defined by the first order systems 

43 ' of differential equations are considered. 

1 Introduction 

The first order polynomial systems of differential equations 

> 

m ; 

depending from the parameters a, b play an important role in various branches of modern mathematics 
and its applications. 
, However even in the case of the system of two equations 

— = k + ax + by + cx 2 + exy + fy 2 , — ^ = I + mx + ny + px 2 + qxy + ry 2 (2) 
as as 

there are many unsolved problems. 

The most famous of them is the 16-th Hilbert problem about the quantity and position of the limit 
cycles in such type of the systems. 

The spatial first order system of differential equations 



=P(?,y,z), ^j- = Q(x,y,z), ^- = R(x,y,z) (3) 
as as as 

with the functions P, Q, R polynomial on variables x, y, z are still more complicated object for the 
studying of their properties. 

As example the studying of simplest spatial systems of equations such as the Lorenz equations 

dx dy dz 

~fo =a v)~ x h -j^=rx-y-xz, — = xy~bz (4) 



or the Rossler system 



dx dy dz 

— = —y — z, —— — x + ay, —— — bx — cz + xz (5) 
ds ds ds 



having chaotic behavior at some values of parameters represent the difficult task. 

The systems of the first order differential equations are not suitable object of consideration from the 
usually point of Riemann geometry. 
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The systems of the second order differential equations in form 

d 2 x l _• . , dx k dx j 

are best suited to do that. 

They can be considered as geodesies of the affinely connected space M in local coordinates x k . The 
values Uj k = U k j are the coefficients of affine connections on M. 

With the help of such coefficients can be constructed curvature tensor and others geometrical objects 
defined on variety M. 

There are many possibilities to present a given system of the first order of equations in the form of 

(6). 

One of them is a following presentation. 
For the system 

£-«••»>. § = m 

after differentiation with respect to parameter s we get the second order system of differential equations 
of the form (6) 

d 2 x 1 dx dy dx 

i 

d 2 y _ 1 (n dx dy dy 

Such type of the system contains the integral curves of the system (7) as part of its solutions and 
can be considered as the equations of geodesies of two dimensional space M 2 (x,y) equipped by affine 
connections with coefficients 

TT 1 — Px TT 1 — Py TT 2 — ® x TT 2 - ® v (Q) 

nn---p, n i2--^ n i2--2g' ^-—q- {9 > 

It is apparent that the properties of the system (7) have an influence on geometry of the variety 
M 2 (x,y). 

Remark that the system (8) is equivalent the second order differential equation 



= M Q /P>,(f) 2 + M W P>^ 



which has the solution in form 

dy _Q 
dx~ P' 

In this connection we can use the second order system of differential equation 

S-^(S)'. S^>.(S) 2 <10 » 

for the studying of the properties of the first order system of equations (7). 

By analogy can be written the spatial system of the first order differential equations. 



2 The Riemann extension of affinely connected space 

For the studying of the geometry of the equations like (8) we apply the notion of the Riemann extension 
of nonricmannian space which was used earlier in [1, 2, 3]. 
Remind basic properties of this construction. 
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With help of the coefficients of affinc connection of a given n-dimensional space can be introduced 
2n-dimensional Riemann space D 2n in local coordinates (x l , having the metric of form 

2n ds 2 = -2U^ j (x l )^ k dx i dx j + 2d^ k dx k (11) 

where ^k are the additional coordinates. 

The important property of such type metric is that the geodesic equations of metric (11) decomposes 
into two parts 

x k + = o, (12) 

and 

1 Jt+lt kji x>x*9i = 0, (13) 

where 

S^k = d^_k u i , T , dxi 

ds ~ ds jk 1 ds 

and R\ji are the curvature tensor of n-dimensional space with a given affine connection. 

The first part (12) of the full system is the system of equations for geodesic of basic space with local 
coordinates x l and it do not contains the supplementary coordinates "J/fc. 

The second part (13) of the system has the form of linear N x N matrix system of second order 
ODE's for supplementary coordinates 

d 2 $ d^ 

-^ + A { s)- + B { s)* = 0. (14) 
Remark that the full system of geodesies has the first integral 

which is equivalent to the relation 

2*fc-j- = ^s + /i (16) 
ds 

where /z, v are parameters. 

It is important to note that the geometry of extended space connects with geometry of basic space. 
For example the property of the space to be Ricci-flat Rij = or symmetrical Rijki-m = keeps also 
for the extended space. 

It is important to note that for extended space having the metric (11) all scalar curvature invariants 
are vanished. 

As consequence the properties of linear system of equation (13-14) depending from the the invariants 
of N x N matrix-function 

e = b- 1 -*±-1a* 

2 ds 4 

under change of the coordinates can be of used for that. 

The first applications the notion of extended spaces for the studying of nonlinear second order 
differential equations connected with nonlinear dynamical systems have been considered in the works 
of author [1, 2, 3]. 



3 The geometry of planar system of equations in form (8) 

We shall consider from geometrical point the properties of planar systems of the equations (7). 

To do this we use the second order system (8) having the solutions of the system (7) as the part of 
their own solutions. 
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Now the system (8) can be considered as the geodesies of two dimensional space with afHne con- 
nections and we extend this space up to the four-dimensional space by the introducing two additional 
coordinates ipi = z, ip2 = t. 

The Riemann metric of the four-dimensional extended space is defined by 



ds A = 2z-0-dx 2 + (2z-± + t ^±)dxdy + 2t-^-dy 2 + 2dxdz + 2dydt. (17) 

Let us consider the basic geometric characteristics of a given metric. 
The nonzero components of the Ricci tensor are 

Rw = 

-3 (&Q{x, y)) 2 P(x, y) + 2 (f?Q(x, y)) P(x, y)Q(x, y) + 2 (JLp(x, y)) Q(x, y)&Q{x, y) 



= 1/2- 



P(x,y) (Q(x,y)f 



- {& p ^,y)) P(x,y) (Q(x,y)f - (&Q(x,y)) (%;P(x,yj) P(x,y)Q(x,y) 
R 12 = -l/2 — ± '- 



(P(x, y )f(Q(x,y)f 



= 1/2 



(mPfry)) (lyP^y)) (Q(x,y)) 2 - (sf^Qfoy)) (P(x, y )f Q(x, y ) 

~ 1/2 (P(x, y )f(Q(x,y)f 

(Q(x,y)f ' 
R22 = 

-3 (j^P{x,y)) Q{x,y) + 2 (jfcP{x,yj)P{x,y)Q{x,y) + 2 (^Q(x,y)) P(x,y)^P(x,y) 

______ . 



As was mentioned above all scalar invariants constructed from the curvature tensor and its covariant 
derivatives 

p = R tJ R^ =0, q = R tjkl R^ kl =0... 

are vanish. 

The full system of geodesies of the metric (17) consists from the two groups of equations 

l_ (£p {x , y) )(*- x{s) ) 2 (%P(x iy ))(£x( S ))£y( S ) 

dsi X[S) P(x,y) P(x,y) 

l_ (&Q(*,g))(£y(«)) a (&Q(s,g))(£s («))&/(*) 

ds 2 [ ' Q(x,y) Q(x,y) 

and 



d: 



((P (x, y)f ( J|Q (x, y)f + Q (x, y) (&Q {x, y)) (jj^Q (x, y)) (P (x, y)) 2 ) t 
^ X °' " Q(x,y) (P(x,y) f 

((P (x, y)f Q (x, y) (£Q (x, y)) J^P (x, y) + 2 (P (x, y)f (Q (x, y)f ^P (x, y)) 



Q(x,y) (P(x,y)Y 

(Q (*> V)f {$P (x, y)) * + (P (x, y)f Q (x, y) £,P (x, y) + (P (x, y)f Q (x, y) (£P (x, y)) 2 ) z 

Q(x,y) (P(x,y)f 



+ - — -2 —+ 



+ ~ — " - 2 — + 



3 THE GEOMETRY OF PLANAR SYSTEM OF EQUATIONS IN FORM (8) 



((Q {x, y)f (x, y)) P (x, y) + (Q (x, y)) 2 (jfcQ (x, yj) P (x, y) J^P (x, y)) z 

Q(x,y) (P(x,y)f 



8 \d (2(£P(x,y))(P(x,y)) 2 Q(x,y)+Up(x,y))(Q(x,y)) 2 P(x,y))£z(s) 



dx^" "7 ds-'"' ' Q(x,y)(P(x,y)) 2 

= 0, 



cP ((P (x, y)f (&P (x, y)) Q (x, y) J^Q (x, y) + (P (x, y)) 3 (£,Q (x, yj) Q (x, yj) t 

I (s) + — - - - h 

ds2 (Q(x,y)) 2 P(x,y) 

(- (P (x, y)f (£Q (x, y)f + (Q (x, y)f P (x, y) J^Q (x, y) + 2 (P (x, y)f (Q (x, y)f (x, yj) 

(Q (x,y)) 2 P(x,y) 

P (x, y) (Q (x, y)f (^P (x, yj) £Q (x, y) + (Q (x, y)f P (x, y) (j^Q (x, y))^ t 

(Q (x,y)) 2 P(x,y) + 
(Q (x, y)f(^P (x, y)Y + P (x, y) (£P (x, y)) (gP (x, yj) (Q (x, y)) 2 ] z 



+ (Q(x,y)) 2 P(x,y) + 

(2 (Q (a;, y)f (j^Q (x, yj) P (x, y) + Q (x, y) (P (x, y)) 2 £Q (x, yj) ±t (s) 
+ (Q (x,y)) 2 P(x,y) + 

+ ^y P ^)Ts Z ^= - 

Last of the two equations looks as the system of 2 x 2 linear matrix equations with the coefficients 
Aij , Bij depending from the solutions of the first system of equations ( this resemble the situation with 
the soliton theory!). 

d 2 z dz , dt _ _ 

— + A n — + A 12 — + B u z + B 12 t = 0, 

ds z ds ds 

d 2 t . dz A dt „ „ 

^ + A 21 —+A 22 — + B 21 z + B 22 t = 0. (18) 

Using the first integral of geodesies of the metric (16) 

dx dy ,„ 
z— + t-f = us 2 + u 
ds ds 

we can get from the system (18) two independent linear second order differential equations for the each 
variable z(s) and t(s) 

d 2 z ., .dz „. , 
- + A(s)- + B(s)z = 0, 

and 

d 2 t .dt „, . 

^ + C(s)- + E(s)t = 0. 

Let us consider some examples. 
1. The system of equations 

dx dy 
ds ^' ds 
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corresponds the second order system of equations 

ds 2 () y{s) ~ U ' 

£ (4-x(a)) 4-y(a) 

as^ x (s) 

From the geodesies of extended space we get the equations for coordinate z 

d 2 (x(,)) 2 z( S ) (2y( S )(x( S )) 2 + 2(y( S )) 3 )fz( S ) 

— rZ 8 ~ 7^ + 1/2-^ s '- hl/° ; 



and 



, z (s) - v K " ; ' + 1/2 = i + 1/2 -f- = 

ds 2 (y( s )f x(s)(y(s)f x(s) 

d 2 , , (y(s)) 2 t(.s) (2(x( S )) 3 + 2x( S )( 2/ ( S )) 2 )^ (s) 

t s - V ^ w; 9 W - 1/2 ^ , '- + 1/2 = 

ds 2 (x(s)) 2 y(s){x(s)) 2 y(s) 

for coordinate t. 

4 The geometry of planar system of equations in form (10) 

The system (10) is more simple object for the studying of the properties of the system (7) from geo- 
metrical point of view. 

Really, the metric of extended space defined by the system (10) can be written in form 

W = _ 2z &^y) dx 2 + 2t d{K !f> y) dy 2 + 2dxdz + 2dydt (19) 
ox ox 

where z,t are the supplementary coordinates and the function K(x,y) — \n(Q/P) is determined from 
the relations 

dlJ , T r, N N 

— = exp{K(x,y)) 

or 

dy , TJ ., ..dx 

- = eMK(x,y))-. 

Contrary to the case (17) the Ricmannian space with the metric (19) is a Ricci-flat 

Rij = 0. 

Its geodesies are defined by the system of equations 

d 2 x dK(x,y) f dx _ . 
ds 2 dx \ ds 

d 2 y dK(x,y) ( dy 



ds 2 dy V^ s / 

5 Rigorous approach at the geometry of the planar systems 

We consider the system of paths of two-dimensional space S 2 m form 

x + YvX^xf + 2Ii\ 2 xy + U\ 2 (y) 2 = 0, 

y + n 2 1 (i) 2 + 2n 2 2 xy + n 2 2 (y) 2 = 0, (20) 
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where coefficients IT* = IT^. 

The Riemann extension of the space S2 is determined by the metric 

4 ds 2 = -2^n} 1 dx 2 - 2m 2 n dx 2 - 4zU\ 2 dxdy - 4tU 2 12 dxdy - 2zli\ 2 dy 2 - 2tTi 2 22 dy 2 . (21) 

A necessary condition for the equations (20) to admit a first integral 

, , dx l 
a i{ x ,y)—r~ = const 
ds 

is 

where 

_ 8ch _ rk 

' j ~ i).r' ij ' 

and r* are the Christoffel symbols of the metric (21) . 

We apply this conditions for determination of the coefficients of equations using the vector in 
form 

a« = [Q(x,y), -P(x,y), 0,0]. 
This means that the first order of equation 

dy = Q(x,y) 
dx P(x, y) 

or 

Q(x,y)dx - P(x,y)dy = 

is an integral of the paths equations. 

Using such condition it is possible to find only three coefficients of affine connections T^-. 
For determination of the rest coefficients we use yet another the first order equation 



(22) 



dy = y(y- i) 

dx x(x — 1) : 
with the first integral 

C(x-l) 

The equation (22) plays an important role in theory of of the planar first order system of equations 
(2) (W). 

In result the coefficients IT^- of the paths equation are defined by 
n i _ (&Q(x,y))x(x-l) 



y 2 P(x, y) - yP(x, y) + Q(x, y)x 2 - Q(x, y)x ' 

n} 2 = 

2 yP(x, y) - 2 P(x, y) + 2 xP(x, y) + (jfcQ{x, y)) x 2 - (jfcQ{x, y)) x - (&P{x, y)) x 2 + (|?(x, y)) x 



= 1/2 

y 2 P{x, y) - yP(x, y) + Q(x, y)x 2 - Q(x, y)x 



(^P(x,y))x(x-l) 



n 1 = 

y 2 P{x, y) - yP{x, y) + Q(x, y)x 2 - Q(x, y)x 
and corresponding expressions for the coefficients nf 1 ,nf 2 ,n| 2 - 
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In result the geodesies take the form 

v 2 



+ 



ds 2 y 2 P(x, y) - yP(x, y) + Q(x, y)x 2 - Q(x, y)x 

(2 yP(x, y) 2 P(x, y) + 2 xP(x, y) + (^q) x 2 - ($q) x- (j|P) x 2 + x) (fajaj) £y(s) 

y 2 P -yP + Qx 2 - Qx 

(*-l)(&(«)) a 
y 2 P - yP + Qx 2 - 

(s) _ (&Q(*,y))y(i/-i)(£*(*)) 2 

ds 2 y 2 P(x, y) - yP(x, y) + Q{x, y)x 2 - Q(x, y)x 

(-(l y Q)y 2 + (^Q)y + 2Qy-2Q + 2Qx + y 2 f x P- y&p) (£x(s)) ±y(s) 

+ yip -yP + Qx 2 -Qx + 

(lp( X ,y))y(y-l)(*-y( S )) 2 

_l i-i L — (23) 

y 2 P{x, y) - yP(x, y) + Q{x, y)x 2 - Q(x, y)x 

+ A (s)f s + B(s)^ s + C(s)z(s) + E(s)t(s) = 0, 

*1 + F(s)^- + H(s)§- + K(s)z(s) + L(s)t(s) = 
as as as 

with the coefficients depending from the parameter s and the functions Q(x,y), P(x,y). 
Remark that last two equations are reduced at the independent equations 

d 2 7 d 7 

- 2+ M [s) - + N {s)z{s) = , 

and 

with the help of the first integral of geodesies 



and 



, . dx , . dy s 
Z ^ + t ^Is- a 2-^° 



of the metric (21). 



6 The second order ODE's cubic on the first derivative in the- 
ory of the planar systems 

The first two equations of geodesic of the metric (23) are equivalent to the one second order differential 
equation 

<tp(x,y))x 2 -(£ y P(x,y))x) (£y{x)f 



y(x) -\ — - - - - - - h 

dx 2 y 2 P(x, y) — yP(x, y) + Q(x, y)x 2 — Q(x, y)x 

((l x P-l v Q)* 2 +{i;Q-l x P-2p) *+(%p) y 2 +{-2P-*-yP)y+2P)(£y(*)) 

y 2 P-yP+Qx 2 -Q(x,y)x 



2 

-+ 
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+ 



y 2 P -yP + Qx 2 

(f x Q(^y))y 2 + (f x Q(^y))y 

y 2 P(x, y) - yP{x, y) + Q(x, y)x 2 - Q{x, y)x 
The equation (24) has the equation 

dy _ Q(x,y) 



+ 77r57Z \ dx i rui x ^\\.2 Z„ — = °- ( 24 ) 



dx P{x, y) 

as particular integral and the function 

>(«) - ^ 05) 

as the first integral. 

In result we get two-parametric solutions of the second order ODE from the one-parametric solutions 
of a given first order equation. 

This fact allow us to describe the properties of the first order system (2) with respect the values of 
parameter C of the solution (25). 

Let us consider some example. 

After substitution of the relation (25) into the equation (24) one gets the expression 

a(x, y)C b + p(x, y)C 4 + j(x, y)C 3 + S(x, y)C 2 + e(x, y)C + n(x, y) = (26) 

, where 

a(x,y) = {a!2 + bl2)y 2 + 
+ (bl + b2 + {2 a22 + 2 all + 2 b22 + 2 bll) x + a2 + al) y+(al2 + bl2) x 2 + (al + b2 + bl +a2)x+ 

+2a0 + 2b0, 



/3(x, y) = 2 b22 y 3 + ((-5 al2 -bl2 -2 b22) x + 2b2- bl2) y 2 + 
+ (H a.22-2, b2-bal-3a2-2 bll -6 b22-bl)x-a2) + 
+ ((-Abll + al2-l0all-bl2-Aa22)x 2 + 2b0-2b2-bl)y+ 
+ (-2 al2 + 2 all)x 3 +{-2 al2-2 bl-2al -3 bl2 -2 a2) .t 2 -4 b0+{-6 aO-2 bl -4 bO-3 b2-al-2 a2) x- 

-2 aO, 



j(x, y) = {-4xb22 - 2 b22) y 3 + ((-bl2 + 2b22 + 10 al2) x 2 - 2 b2 + (bl2 + 4 b22 - 4 b2) x) y 2 + 
+ (b2 + (-3 al2 + 2 a22+2 bll + bl2 + 20 all)x 3 + {<6 b22 + 2 a22 + bl +3 a2-A bO + 6 b2) x) y+ 
+ ((-bl +10 al+3 a2 + 2 b2 + A bll +2 bl2-al2 + 8 a22) x 2 -2 bO) y+ 
+ (al2 -6all)x 4 + (4al2 + a2 -2 all + bl) x 3 + (bl + 6 aO + 3 b2 + a2 + 8 bO) x+ 
+ (al2 +6a0 + Abl + 2 al +2b0 + 3bl2 + Aa2)x 2 + 2b0, 



5(x,y) = (2x 2 b22 +Axb22)y 3 + ({bl2 - 10 al2)x 3 + (2b2-4b22 + bl2)x 2 + (4 b2 -2b22)x) y 2 + 

+ ((bl -Aa22 + 3al2 -2 bll - a2 - 10 al - 2 bl2) x 3 + (-bl2 - Ab2 - Aa22 - 3 a2 + bl +2b0) x 2 ) y+ 



6 THE SECOND ORDER ODE 'S CUBIC ON THE FIRST DERIVATIVE IN THEORY OF THE PLANAR SYSTEMSW 



+ ((3 al2 - 20 all ) x 4 + (-3 b2 + 4 60 - 2 b22) x) y+ 
+6x 5 all + (6 all + 2 al -2al2)x 4 + (-4 60 - bl2 -2 61 - 2 a2 - 6 aO) x 2 + 
+ (-2al2 -2a0 -2a2 -2bl)x 3 + {-AbO - b2)x, 



e(x,y) = -2y 3 x 2 b22 + (-x 3 b!2 + hx 4 al2 + (-2 b2 + 2 b22) x 2 ) y 2 + 
+ ((5 al -3 al2) x 4 + (10 all - al2) x 5 + {-bl + bl2 + a2 + 2 a22) x 3 + {2b2-2 bO) x 2 ) y-2 x 6 all + 
+ {-al - 6 all ) x 5 + (bl + 2 aO + a2) x 3 + {a!2 -2al)x 4 + 2 x 2 b0, 



fi(x,y) = -y 2 x 5 al2 + (-2x 6 all + (a!2 - al) x 5 ) y + 2 x 6 all + x 5 al . 
The substitution of the expression (25) into the relation (26) give us the expression 

A(x)C 6 + B(x)C 5 + E{x)C 4 + F(x)C 3 + H(x)C 2 + K{x)C + L(x) = 0, (27) 

where 

A(x) = (2 bl2 - 2 b22 - 2 bll + 2 a!2 - 2 a22 - 2 all ) x 2 + 
+ (2 b22 - 2 bl2 - 2 al2 + 2 all +2a22 + 2 bll ) x+ 
+ a2 + al2 + bl +bl2 + b2 + 2a0 + al + 2 bO, 

B(x) = (-8 al2 + 12 all + 4 bll - 4 b22 + 4 a22) x 3 + 
+ (-2 bl +7 al2 + 16 b22 + 2 al + 4 b2 - 2 bll - 10 all - 5 bl2) x 2 + 
+ (-2bl -8b2 -8a0 -4a22 -Aa2 + 2bl2 - 4 a!2 - 14 b22 - 8 bO - 6 al - 2 bll ) x- 

-2 60 -2a0 - bl - a2 - bl2 + 2 b22 

E(x) = (-28 all - 2 bll - 2 a22 - 2 b!2 + 12 al2) x 4 + 
+ (-8 al + 18 all + 6 b22 + Abl - 2 bll + Abl2 - 6 a22 -8al2-A b2) x 3 + 
+ (4 bl + 12 aO + 6 a2 + 13 al + b2 + 6 al2 + 6 a22 + Abll + 4 bl2 - 20 b22 + 10 bO) x 2 + 
+ (4 a2 + 2 bl + 10 bO + 10 b2 + 2 a22 + 16 b22 + 8 aO) x - b2 - 2 b22, 



F(x) = (32 al 1 - 8 al2) x 5 + (-2 bl + 4 a22 + 12 + 3 bl2 + 2 bll + 2 al2 -12 all) x 4 + 

+ (-4 a2 - 6 bl - 12 al - 4 a!2 - 6 bl2 - 8 aO - 2 &£g - 2 611 - 4 60 + 6 6^) x 3 + 
+ (-2 61 - 8 b2 - 14 60 - 12 aO - 61^ - 6 a2 - 4 a22 + 6 6^) x 2 + (-4 b22-2b0-2 b2) x, 



H(x) = (2al2 - 18 all) x 6 + (2al2 -2 all -8al)x 5 + 
+ (-2 a22 + 2a0 - bl2 + a2 + 3 61 + 3 al + al2) x 4 + 
+ (2 61 + 2 612 + 6 60 + 4 a2 - 2 62 + 8 aO + 2 a22) x 3 + (3 b2 + 4 60) a; 2 , 
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K{x) = 4x 7 all + (Sail + 2 al - al2)x 6 + 2x 5 al + (-2 aO - a2 - bl)x i -2x 3 b0, 

L{x) = -2x 7 all -x 6 al. 
The substitutions x = and x = 1 into the (27) lead to the conditions on the value C 

(a2 + al2 + bl + bl2 + b2 + 2 aO + al + 2 bO) C 2 + 

+ (-2 bO -2a0 - bl - a2 - bl2 + 2 b22) C - b2 -2b22 = 0, (28) 

(a2 + al2 + bl + bl2 + b2 + 2 aO + al + 2 bO) C 2 + 
+ (2 all - 2 bO - bl - 2 aO - al2 - a2)C -2 all - al = 0. (29) 
The substitution (x = C, y = 1 — C) into the (26) lead to the conditions on the value C 

(bl2-2b22)C + 2b22 + b2 = 0. (30) 

After substitution (y = 1 — x, x = 1, C = 1/C1) we get 

{al + 2 all) CI 2 + (a!2 - 2 all + bl + a2 + 2a0 + 2b0) CI - al - a2 - al2 - 2 bO- 

-bl2 -bl -b2 -2a0 = 0, (31) 
and the substitution (y = 1 — x, x = 0, C = I /CI) lead to the condition 

(-2 b22 - b2) CI 2 + (-2 bO - bl2 - bl - a2 - 2 aO + 2 b22) CI + al + a2 + al2 + 2 b0+ 

+ bl2 + bl + b2 + 2a0 = 0. (32) 

Remark that for the comparison with results of the article ([4]) in all above formulas the equation 
(2) was presented in the form 

dy aO + alx + a2y + allx 2 + al2xy + a22y 2 



dx bO + blx + b2y + Mix 2 + bl2xy + b22y 2 ' 

Remark 1 In the famous article ([4]) was developed the approach to the studying of the problem of the 
limit cycles of the equations (2). 

Let us remind the basic facts of the Petrovsky-Landis theory. 

For every closed curve of the system (2) the solution 

V{ >~ (x-C) 

of equation 

dy = y{y-l) 
dx x(x — 1) 

is corresponded. 

At the same time the value C satisfies the algebraic equations 



5>„(^)C" = 0, 
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where the coefficients a n ((ii) are dependent from the parameters of equation (2). 
This equation arises from the condition 

f (x - C) 2 [x(x - l)Q(x, y) + y(y - lJPjx, y)} ^ = Q 

Jc X 3 (X-1) 3 

The substitution of the function y from the (33) into this expression and calculation of the residues 
with regard the points 0, 1, C lead to the equations on the parameter C. 

According the ([4]) general quantity of the values C defined by such equations is equal 14 and this 
number coincide with the quantity of closed solutions determined by the equations (2). 

As it was shown in ([4]) 11 curves from 14 can be transformed into the vicinity of the particular 
points of the equation (30). 

As result only three closed curves do not be transformed into the vicinity of the particular points and 
the quantity of the limit cycles defined by the equation (2) was found equal three. 

In spite of the fact that the statement of the article ([4]) about the quantity of the limit cycles defined 
by the system (2) was fallacious its approach is useful from point of our consideration. 

In fact the conditions on the value C given by the (28)- (32) are the same with those which was used 
in the article ([4]) for estimation of the quantity of the limit cycles in quadratic polynomial system. 

It can be shown that all conditions on residues of functions considered in the article ([4]) are followed 
from the (28)- (30). 



7 The examples 

We apply the second order ODE (24) for the studying of the properties of the first order ODE's 

dy = Q(x,y) 
dx P{x, y) 

1. Let us consider the system 

y = 8 — 3a — lAax — 2axy — 8y 2 , 

x = 2 + 4x- 4ax 2 + 12xy. (34) 
Corresponding equation (24) for this system becomes 

d 2 

(-12y(a;) 3 a;+(-2+4aa; 2 + 8.T 2 ) y(x) 2 + (2ax 3 + 2-6ax 2 +4x) y(x) -3 ax-8 x 2 + 8 x- 11 ax 2 + U ax 3 ) — y{ 

+ (-Ay(x)x 2 -4 + 4 x 2 + 2 ax 2 + 24y(x)x - 8 ax 2 y(x) -2 ax 3 + 12 (y(x) f x + Ay(x)^j (j^v( x ) 

+ (-12y(x) 3 + (10ax + 16x + 8)y(x) 2 + (6a + 20ax + 2ax 2 -l2) y(x) + l6-8ax + Uax 2 -6a-l6x) ^-y{x)+ 

dx 

+14 ay(x) - 12 a (y(x)) 2 - 2 a (y{x)) 3 = (35) 

and the function 

dy(x) 8 — 3a— 14ax — 2axy — 8y 2 
dx ~~ 2 + 4x- 4ax 2 + I2xy 

is solution of this equation. 
Moreover the function 

dy{x) = y(y- 1) 
dx x(x — 1) 
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also is solution of the equation (35). 

Now we shall get another particular solutions of our equation (35). 
The substitution 

dy = M(x,y) 
dx N(x, y) 

into the equation (35) give us the relation between the function M(x,y), N(x,y). 

In spite of the fact that such relation is very complicated one can get with help of the MAPLE-6 
some solutions of them. 

For example, the presentation of the functions M(x,y), N(x,y) in form 

M(x, y) = A(x) + B(x)y + C(x)y 2 , N(x, y) = E(x) + F(x)y + H(x)y 2 , 

where 

A(x) = 2x- 3ax 2 - 1, B(x) = -1, C(x) = -2x, 
E(x)=x, F(x) = 2x 2 , H(x)=0, 

corresponds the system 

dy(x) 2x — 3ax 2 — 1 — y — 2xy 2 
dx x + 2x 2 y 

compatible with the equation 

d 8 — 3 a — 14 ax — 2axy(x) — 8 y{x) 2 

dx V ^ X ' ~ 2 + 4 x - 4 ax 2 + 12 yjxjx ' 

In result we get the solution 

1/4 + x - x 2 + ax 3 + xy{x) + x 2 y{x) 2 = 

which present the limit cycle at the some value of the parameter a. 

2. Let us consider the system 

y = x + 2y + Axy + (2 + 3a) y 2 , 

x = 5x + 6x 2 + 4(1 + a)xy + ay 2 . (37) 
Corresponding equation (24) for this system becomes 

(2x 2 ay — Ax 2 a — 2xay + 4x 3 a + 4x 3 — 4x 2 ) (-j^y(x) 
A(x,y) 



, 3 



, 2 



(-6xy - 20yx 2 - 6xy 2 a - 3i 2 - lQx 2 ay - Ay 2 x - Ax 3 + 7x + 6xay) {-^y{x))' 
+ A(x,y) + 

(6 xy 2 a - 4 ay 2 - 7 y + 3 y 2 + 20 y 2 x - x + x 2 + 4 y 3 - 6 xy + 4 yx 2 + 4 ay 3 ) j^y{x) 
+ A(a;,y) + 

d 2 y + 3y 2 -4y 3 

+ dx^ yix) + A(x,y) = °' (38) 

where 

A(x,y) = y 4 a+(-a+4:X + 4:xa) y 3 +(-7 xa-x + 3x 2 a+8x 2 ) y 2 +(-8x 2 -7 x+Ax 3 ) y-x 2 + x 3 . 
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The first order equation 

dyjx) = x + 2y + 4xy+{2 + 3a)y 2 
dx 5x + 6x 2 + 4(1 + a)xy + ay 2 ' 

is particular integral of the second order equation (38). 
The function y(x) defined by the relation 

x 2 + x 3 + x 2 y(x) + 2axy(x) 2 + 2axy(x) 3 + a 2 y{x) A = (40) 

satisfies both equations (38) and (39). 

Beyond this point the first order equation 

dy{x) _ 2x + 3x 2 + 2xy + 2ay 2 + 2ay 3 
dx x 2 + 4axy + 6axy 2 + Aa 2 y 3 

resulting from the (40) by differentiation also satisfies the equations (38). 

The presence of the two particular first order integrals for a given second order differential equation 
can be used for finding one them in evident form. 

In concerned case we use the integral (39) as known and seek for the expression 

dy(x) = M(x,y) 
dx N(x,y 

as the second integral with some functions M(x, y) and N(x, y) ( for example as polynomial on variable 
y with the coefficients depending from the variable x). 

After substitution both particular integrals into a given second order ODE (38) the functions M(x, y) 
and N (x, y) can be find in evident form from the corresponding condition of compatibility. 



3. The next example is the system ([6]) 

y = 15(1 + a)y + 3a(l + a)x 2 - 2(9 + ba)xy + I6y 2 , 

x = 6(l + a)x + 2y~6(2 + a)x 2 + l2xy. (41) 
Corresponding equation (24) for this system becomes 

(-2x+12x 3 -10^ 2 ) (£y(x)) 3 
A(x, y) 

(I%xy-2$yx 2 + l%x 3 -mx 2 -2y 2 -12y 2 x+2lx-l2yxa+l()x 3 a-31ax 2 + l2yx 2 a+21xa+2y) (£y{x)) 2 | 
+ A(x,y) h 

| (6yx 2 a 2 -18 yx 2 -21 y +21 y 2 a-10 y 2 x+l2 y 3 -21 ya-22y 2 xa+%y 2 -4yx 2 a+Mxy + 42yxa) -£y(x) | 

H M^y) 2 h 

d 2 , N 18y 3 -18y 2 + 6xa 2 y-10y 2 a-6xa 2 y 2 + 10y 3 a-6y 2 xa+6yxa , , 

^ y(x)+ - — ■ ■ — wj) ' — ' = °' (42) 

where 

A(a;, y) = (2 + 12 x) y 3 + (-22 x + 4 x 2 - 2 - 6 ax 2 + 6 xa) y 2 + 
+ (-10 x 3 a - 18 x 3 + 31 ax 2 - 21 xa + 45 x 2 - 21 x) y - 3 x 3 a - 3 x 3 a 2 + 3 x 4 a + 3 x 4 a 2 . 

8 The spatial first order system of equations 

Generalization of the Petrovsky-Landis approach to the spatial first order systems of equations may be 
realized by the following way. 
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8.1 

Instead of the spatial first order system of equations 

dx dy dz 

— =P{x,y,z), — = Q{x,y,z), —= R(x,y,z) (43) 

we consider the Pfaff equation 

P(x, y, z)dx + Q(x, y, z)dy + R(x, y, z)dz = (44) 

for the orthogonal trajectories of the system (43). 

The relation (44) can be considered as the linear first integral of geodesies of the the six-dimensional 
space in local coordinates (x, y, z, U, V, W) with the Riemann metric of the form 

6 ds 2 = -2{T\ X U + T 2 n V + T 3 n W)dx 2 - 2(Tl 2 U + T 2 22 V + T\ 2 W)dy 2 - 

-2(Tl 3 U + T 2 33 V + T\ 3 W)dz 2 - A{T\ 2 U + T\ 2 V + T\ 2 W)dxdy- 
-A{T\ 3 U + T\ 3 V + T\ 3 W)dxdz - A{T\ 3 U + T 2 23 V + Tl 3 W)dydz+ 

+2dxdu + 2dydV + 2dzdW, (45) 

where = T^(x,y,z). 

One part of geodesies of the metric (45) has the form of deodesics of the three-dimensional space 
with the affinc connection = F^(x, y, z) 

n x ^ (dx\ 2 j (dy\ 2 , (dz\ 2 , dxdy , dx dz , dy dz 



^ +r "UJ +T22 {fs) +r -UJ +2ri2 ^ +2rl3 ^ +2r ^ =0 ' 

d 2 y ^ 9 / dx\ 2 o / dy\ 2 ~ / dz\ 2 dxdy „ ~ dx dz dy dz 

^ +r "U) +T22 {I) +r -U) +2r ^i +2r "^+ 2r -i^ = °> 

d2z , r.3 ( dx Y , t.3 ( d V\ 2 , r 3 ( dz Y , ot-3 dxd V , ot-3 dx dz , OT -3 ^ dz 



, .-rfihr + r 22 -f + r 33 -r + 2r i 2 ^-f + 2^3— — + 2r^-^— = 

ds 2 V rf -V \ ds J \ds J 12 ds ds u dsds 26 ds ds 

and another part of geodesies is the linear system of the second order differential equations for the 
coordinates U, V, W. 



Remind that if the relation 



is the linear integral of geodesies 



then the conditions 



dx 2 „ 
di— = 
ds 



d 2 x l : dx j dx k _ 
~ds T + ]k mis~ = 

+ dj-i = 



must be satisfy. 

In the space with the metric (45) these conditions allow us to determine only six coefficients 
T k j(x,y, z) from the eighteen one. 

For determination of others coefficients T^{x,y,z) we use the auxiliary equations same with the 
Petrovsky-Landis theory 

y(y — l)dx + x(x — l)dy = 0, z(z — l)dx + x(x — l)dz = 0. 

Every of these relations can be considered as the linear first integral of geodesies of the metric (45) 
and this allow us to determined all coefficients of connection T^(x, y, z). 
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Some of them are 



= 1/2 













(&P{x,y,z))x{x-l) 






P(x,y, z)x 2 - 


P(x, y, z)x - Q(x, y, z)y 2 + Q(x, y, z)y - 
r n (x,y,z) = 


- R(x,y, z)z z + R(x,y, 


z)z 




y(y-l)3-P( X ,y,z) 






P(x,y,z)x 2 - 


- P(x, y, z)x - Q(x, y, z)y 2 + Q(x, y, z)y 
T 3 n (x,y,z) = 


- R(x,y,z)z 2 + R(x,y. 


,z)z 




(&P(x,y,z))z(z-l) 






P(x,y,z)x 2 - 


- P(x, y, z)x - Q(x, y, z)y 2 + Q(x, y, z)y 
F 2 13 (x,y,z) = 


- R(x,y,z)z 2 + R(x,y. 


,z)z 




-^R{x, y, z) (x - x 2 )+2R{x, y,z)(x-l- 


^z)-JLp(x,y,z)(x 2 - 


x)) 



rL(x,y,z) 
^2(x,y,z) 

Tl 2 (x,y,z) 



(P(x, y, z)x 2 - P(x, y, z)x - Q(x, y, z)y 2 + Q(x, y, z)y - R(x, y, z)z 2 + R(x, y, z)z) x{x-\) 

(j^Q{x,y,z)} x(x - 1) 



P(x, y, z)x 2 - P(x, y, z)x - Q{x, y, z)y 2 + Q{x, y, z)y - R{x, y, z)z 2 + R(x, y, z)z ' 

, y(y-i)%Q(x,y,z) 

P(x, y, z)x 2 - P(x, y, z)x - Q(x, y, z)y 2 + Q(x, y, z)y - R(x, y, z)z 2 + R(x, y, z)z ' 

z ( z - l )-§iQ{x,y,z) 

P(x, y, z)x 2 - P(x, y, z)x - Q(x, y, z)y 2 + Q(x, y, z)y - R(x, y, z)z 2 + R(x, y, z)z ' 
x{x-l) (^R{x, y, z) + £Q(x, y, «)) 



Tl 3 (x,y,z) = 1/2 



= -1/2 



P(x, y, z)x 2 - P(x, y, z)x - Q(x, y, z)y 2 + Q(x, y, z)y - R(x, y, z)z 2 + R(x, y, z)z ' 

Th(x, y,z) = 



^R(x,y,z) + f z Q(x,y,z))y(y-l) 



P(x, y, z)x 2 - P(x, y, z)x - Q(x, y, z)y 2 + Q(x, y, z)y - R(x, y, z)z 2 + R{x, y, z)z ' 

^2^{x,y,z) = 

z ((j^R(x,y,z)j z - -^R(x,y,z) + (J^Q(x,y,z)) z - -^Q(x,y,z)j 
1//2 P(x, y, z)x 2 - P(x, y, z)x - Q(x, y, z)y 2 + Q(x, y, z)y - R(x, y, z)z 2 + R(x, y, z)z ' 



T 33(x,y,z) 
?33(x,y,z) = 



r 33(x, y,z) 



x ( x ^ x ) -§i R {x,y,z) 



P(x,y 1 z)x 2 - 


- P(x,y, z)x - 


- Q(x,y, z)y 2 + Q(x,y, z)y - 
(f z R(x,y,z))y(y-l) 


- R(x 1 y,z)z 2 - 


s rR{x,y,z)z' 


P(x,y,z)x 2 


- P{x, y, z)x 


- Q(x,y,z)y 2 + Q(x,y,z)y 
z(z-l)&R(x,y,z) 


- R(x,y,z)z 2 


+ R(x,y,z)z 


P(x,y,z)x 2 


- P(x,y,z)x 


- Q(x,y,z)y 2 + Q(x,y,z)y 


- R(x,y,z)z 2 


+ R(x,y,z)z 



So the six-dimensional space with the metric (45) is suitable geometric object for the studying of 
the properties of the spatial first order system of equations (43). 
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8.2 

Another approach for the studying of the spatial system of equations is connected with consideration 
of the system (43) 

Q(x,y, z)dx — P(x,y, z)dy = 0, R(x,y, z)dx — P(x,y, z)dz = (46) 
together with the relation 

xdx + ydy + zdz = 0, x 2 + y 2 + z 2 = C (47) 

( or another ones) as the linear integrals of geodesies of the six-dimensional space. 

From these conditions the eighteen coefficients of connections r*- can be uniquely determined. 
In fact, we get the following equations 

+ y, *)) g + gfo y, *) + y&gfo y, *)) (jA^f { 

ds 2 yQ{x, y, z) + R(x, y, z)z + P(x, y, z)x 



+ 



0, 



(- (&Pfa y, g) V + (jyQi^ y, zj) y + (l y R{x, y, zj) z) ( 4- s y{s) 
yQ{x, y, z) + R(x, y, z)z + P(x, y, z)x 

((jLp(x,y,z))y-P(x,y,z)) (&y( s )) 2 

yQ{x, y, z) + R(x, y, z)z + P(x, y, z)x 

(j£P(x, y, z)) z - (JLR(x, y, zj) z - (&Q(x, y, z)) y) (£x(s)) £z(s) 
yQ{x, y, z) + R(x, y, z)z + P(x, y, z)x 

({^Pjx, y, zj) z + ylPjx, y, zj) jzjs) ((^Pfr, y, zj) y P(x, y, zj) (fajstf 

yQ(x, y, z) + R(x, y, z)z + P(x, y, z)x yQ(x, y, z) + R(x, y, z)z + P(x, y, z)x 

and analogous equations for the 

5? y(«) + ... = o 

and 

5? ,(-) + ... = o 

In result the six-dimensional space with the geodesies, having the linear integral in form (46)-(47) 
has been constructed. 

9 Projectivization of the planar system 

We apply the geometric approach at the studying of the Pfaff equations (44) connected with a planar 
systems of equations 

g=p(*,V), §=«(*,*) (48) 

after their projectivization. 

In fact after the continuation of the system (46) on the projective plane we get the Pfaff equation 

-zQ(x, y, z)dx + zP(x, y, z)dy + (xQ(x, y, z) - yP(x, y, z))dz = 

where P(x, y, z), Q(x, y, z) are the homogeneous functions constructed with the help of given functions 
p(x,y),q(x,y). 
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As example we consider the equation 

^ = Ax - y - 10x 2 + (5 + 5)xy + y 2 , ^ = x + x 2 + (e - 25)xy 

having at least a four limit cycles at the some restriction on the parameters A, S, e. 
After its projectivization we get the Pfaff equation 

(—y 2 x5 — yXxz + zy 2 — 5 y 2 x + zx 2 + x 3 + x 2 ye — 15 x 2 y — y 3 ) dz+ 

+ (z 2 A x — z 2 y + 5 zxy + zxyS + zy 2 — 10 zx 2 ) dy + (~zxye + 25 zxy — z 2 x — zx 2 ) dx = 0. 

For determination of the connections coefficients of corresponding six-dimensional space we add at 
this equation the relations 

y(y — \)dx + x(x — l)dy = 0, z(z — l)dx + x(x — l)dz = 

and consider all of them as the linear integrals of geodesies. 

In result we get a six-dimensional space whose geodesies have at least a four limit cycles. 

9.1 The Lorenz system 

The next example is the system of equations 

„, . dy . dz 

— =P(x,y,z), — =Q{x,y,z), — = R{x,y,z) 

with some functions P,Q,R. 
We consider the relations 

Q(x, y, z)dx — P(x, y, z)dy = 0, R(x, y, z)dx — P(x, y, z)dz = 

as the linear first integrals of geodesies of corresponding six-dimensional space. 
For determination of the connection coefficients we use the invariant conditions 

rli + r 2 2 + r 3 3 - o, v\ 2 + r 2 2 + r 3 3 = o, r{ 3 + r 2 3 + rj 3 = o, 

and the freedom in the choice of coordinate systems. 

We use a following normalization of coordinate system 

rli = o, r} 2 = — r 23 , r 22 = o, r 2 2 = — r 3 3 , r 33 = o r} 3 = — r 23 . 

In result all connection coefficients of the space can be determined uniquely and we get the expres- 
sions for the connection coefficients 

-Q-g-P + RjLp + P-g-R - 2 P-g-P + P#Q 
T 3 23 (x,y,z) = -1/4-^ ^ dz 



PQ 

r2 i \ Q( x ,y,z)-^P(x,y,z) 



P(x,y,z)R(x,y,z) 

-RQ^-P + (R) 2 #P + PR^-R - 2 QP^-R - 2 PR-g-P + PR-g-Q 
T\ 2 {x,y,z) = 1/4 -^ — ^ <±_dy 



Q(py 

r2f , 2 Qp£p QPf y Q (Of lyP + QRI.p QP&R + 2 (&Q) RP 

1 13\ x , Vi z ) — ^V 4 (P) 2 "/? ' 
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■g-P(x,y,z) 

3 R(x,y,z)£jP(x,y,z) 
P(x,y,z)Q(x,y,z) 

_2 P-S-P - 

T\ 3 {x,y,z) = -l/A 



2 PJLp + Pf y Q + Qf y P + Rf z P + PJLR 



RQ 



m P{x,y,z) 2 g^Q(x,y,z) 



Tl 3 (x,y,z) = - dz " T 2 u (x,y,z) 



R(x,y,z) ' ' ' ' P(x,y,z) 

T 3 n (x,y,z) = a * V ' 



P(x,y,z) 



2 P#- P - P#- Q - Q-S-P + P#-P - P#-P 

r:i,(x.. (A2 ) = i/4. " ' " 



PR 

-PttQ - QttP + R-n-P + P^-R 



rf 3 (a;,y,z) = -1/4 ^ 

The corresponding metric is 

6 ds 2 = -2 (^(ar, y, z)V + T 3 n (x, y, z)W) dx 2 - 

-4 (-T 3 23 {x, y, z)U - T 3 3 (x, y, z)V + T 3 2 (x, y, z)W) dxdy- 
-4 (-T 2 23 (x, y, z)U + T 2 3 (x, y, z)V + T 3 3 (x, y, z)W) dxdz- 

-2 (Tl 2 (x, y, z)U + T 3 2 (x, y, z)W) dy 2 - 
-4 (T^ix, y, z)U + T 2 23 (x, y, z)V + T 3 23 (x, y, z)W) dydz- 
-2 (Tl 3 (x, y, z)U + T 2 33 (x, y, z)V) dz 2 + 2dxdU + 2dydV + 2dzdW. 
Remark that the conditions 

ai-j-,k + RTi^m = 

for the vectors 

a t = [P(x, y, z), -Q(x, y, z), 0, 0, 0, 0] 

and 

a, = [P(x, y, z),0, -R(x, y, z), 0, 0, 0] 

are obeyed for the metric (49). 

In particular case of the Lorenz system of equations 

dx . , dy dz 

— =<j{y-x), —=rx-y-xz, —=xy-bz, 

ds ds ds 

we have 

P(x,y,z) = a(y - x), Q(x,y,z) = rx - y - xz, R(x,y, z) = xy - bz, 
where cr, r, b are the parameters and geodesies of the metric (49) are equivalent to the expressions 

*_ y ( x ) ° iM^)l + i/ 2 (x + xb-2x* + rx-zx + 2<Ty-yb-2y)* {±y{x)f &z(x) 
dx 2 —rx + y + zx rx 2 y — yzx 2 — y 2 x + xbz 2 — zxrb + ybz 

(x + xb - 2 xa + rx - zx + 2 a y - yb - 2 y) (£y(x)) £z(x) 
bzx — ybz — yx 2 + y 2 x 
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-1/2 
-1/2 



(— yb + xb + zx — 2 xo ■ 



2" V) (&(*))' 



rxy — y 2 — yzx — rx 2 + xy + zx 2 

(-yb + xb + zx - x - rx + 2y) ^y(x) z-r 

(y 2 - 2 xy + x 2 ) a a (-y + x) 



+ 



^2 + 1/2 (x + xb-2xn + rx- zx + 2ay-yb-2y)n (£y(x)) (£z(x)) 2 
dx 2 rx 2 y — yzx 2 — y 2 x + xbz 2 — zxrb + ybz 

°_ (£y( x )) 2 S z( +) _ (-V b + xb + zx-2x<7 + x-rx + 2ay) (£y(x)) £z(x) ^ 
—rx + y + zx rxy — y 2 — yzx — rx 2 + xy + zx 2 



(-xy + bz) (sj/M) 1 , 2 (x + xb - 2 xcr + rx - zx + 2 a y - yb - 2 y) (-^z{x))' 

rxy — y 2 — yzx — rx 2 + xy + zx 2 bzx — ybz — yx 2 + y 2 x 

(xbz 2 + ((-3x 2 -b 2 + 2ab) y + 2 x 3 + (b 2 -2 a b+b -rb) x) z + (-2a + b-2)xy 2 + (3r-b+2a + l)x 2 y-2x 3 r) £y(x) 



+1/2 

(— y + x) (—rx + y + zx)a 



(—yb + xb + zx — x — rx + 2 y) -^z(x) y 

(y 2 — 2 xy + x 2 ) a (—y + x)a 



and the second order linear system differential equations for the coordinates U(s),V(s), W(s). 
These equations have the solutions in form of the first order differential equations 

dy rx — y — xz dz xy — bz 



dx a(y — x) dx u(y — x) ' 

Geometrical properties of the metric (49) are dependent from the parameters a,b,r. 
For example the component R zz the Ricci tensor Rij look as 

1 x (-zx +(2a-b-2)y+{l + b-2a + r)x) 
(—xy + bz) (—y + x) (—rx + y + zx) 

and a more complicated expressions for another components. 

Finally we present the expression for the density of the Chern-Simons invariant of affine connection 
for the space defined by the Lorenz system of equation. 

In case of three dimensional space it is defined as ([6]) 

cs(Y) = e« fc (r? + 2 -v% r« r r fep ) (50) 

where e yfe is a Levi-Civita symbols. 

For the six-dimensional metric we get the expression 

4cr (—y + x) A (—rx + y + zx) 2 (—xy + bz) 2 CS(T) = 

= -9z 5 x 4 b+((l0x 5 + (lib 2 + 48 ab- 40 b) x 3 ) y + (-11 b 2 + 4 b + 36 rb - 25 cr b) x 4 - 23 ay 2 x 2 b - x e ) z 4 + 

+M3(x, y)z 3 + MA(x, y)z 2 + M5(x, y)z + M6(x, y), (51) 

where 

M3(x, y) = (4 r - 2 a) x 6 + (18 b - 40 r - 8 + 32 a) yx 5 + 
+ ((-18 b + 44 - 56 cr) y 2 - 14 a b 2 - 54 br 2 + 12 b 2 + 7 b - 18 a 2 b - 4 a b + 75 a rb + 5 b 3 - 12 rb + 26 rb 2 ) x 4 + 
+ (26 ay 3 + (-144 orb - 2 b + 80 a 2 b + 120 rb - 66 a b - 10 b 3 -50b 2 - 26rb 2 + Uab 2 )y) x 3 + 
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+ (-106cr 2 6+ Uab 2 + 5b 3 + 38b 2 + 69 a rb- 59 6+ 138crfe) y 2 x 2 + (Ua 2 b-68crb- Uab 2 )y 3 x. 

The functions M4(x,y) — M6(x,y) are dependent from the parameters and are too cumbersome. 
Remark that in the case x = y for the right part of the (51) we get the expression 

-9y 4 (z + l-r) 4 (bz-y 2 ). 

which is equal to zero on the stationary points of the Lorenz model 

z = r - 1, y = y/b(r - 1). 

More detail information on the properties of the six-dimensional Riemann space connected with 
the Lorenz system of equations can be obtained with the help of the studying of the linear system of 
equations for the coordinates U, V, W. 

The studying of the equation 

g l > ViVj A k - R l k A x + XA k = 
for the eigenvalues A of the dc Rham operator, defined on the 1-forms 

A(x,y,z) = A l (x,y,z)dx t 

of manifold can be also useful for that. 
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